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. Abstract. We investigate, in a set of 3D numerical simulations of driven, 

Q ' magnetized, isothermal, and self-gravitating molecular clouds (MCs), the statis- 

D , tical correlations between the energy ratios (thermal / gravity, and kinetic / gravity) 

^"*^ ■ of clumps and cores (CCs) identified in the simulations and gravitational bind- 

ing indicators commonly used in observational studies such as the Jeans number, 
J c , and the virial parameter, a V i r . In the energy ratios, we consider the surface 
energy terms which account for the effects of the environment on the clump 
gravitational boundedness. We find that: a) J c and the thermal/gravitational 
^ , energy ratios are well correlated, b) a V i r and the (thermal+kinetic) /gravity or 

iy-j 1 kinetic/gravity energy ratios are poorly correlated, additionally affected by the 

ambiguity of the compressive or dispersive effect of the velocity field. This result 
\Q , suggest that the use of a V i r estimates in the observations is only useful to assess 

CN) ' the kinetic+thermal energy content of a CC and not its gravitational bounded- 

ness. Finally, we discuss briefly the possibility of measuring the kinetic energy 
. surface term directly in the observations. 
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1. Introduction 

One of the most challenging issues in the field of star formation is to assess the 
gravitational boundedness of clumps and cores (CCs) which are prone to star 
formation. Constraining the masses and scales on which CCs are gravitationally 
bound sets important constraints on the star formation efficiency in molecular 
clouds (MCs). Two quantities commonly used in the observations in order to 
assess the gravitational boundedness of CCs are the Jeans number J c and the 
virial parameter a v i r . J c measures the importance of the thermal support against 
gravity and is defined as: J c = R c /Lj )C , where L j jC = (Ttc^/Gpc) 1 / 2 is the 
mean Jeans length in the core, p c the average density, and c s the thermal sound 
speed, and R c the characteristic size. The virial parameter measures the ratio 
of thermal+kinetic to gravitational energies (e.g., Bertoldi & McKee 1992) and 
is defined ELS. Q'uiv 

= (5 a 2 R c /G M c ), where M c and a are the mass and one 
dimensional density-weighted velocity dispersion of the CC. 

Another method for evaluating the gravitational boundedness of CCs is by 
using the virial theorem (VT). As shown in Dib et al. (2006a), the surface 
energy terms in the VT are equally important for the overall virial balance. 
Thus, in order to compare the gravitational term W in the VT to the net ef- 
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Figure 1. (Left) Relationship between J c and the ratio of thermal to grav- 
itational energies for CCs in a strongly supercritical cloud model. All clumps 
have (E t h — Wi) > 0. The full line is a fit to the data when the thermal 
surface energy terms are taken into account (stars) and the dot-dashed line 
when they are not (triangles). (Right) Relationship between a V i r and the 
ratio of the kinetic+thermal to gravitational energies in the same model. 
CCs are cataloged by whether kinetic+thermal energies provide a net sup- 
port to the clump (Ek + E t h — t& — Tth > 0, stars), or a net compressive effect 
(Ek + E t h — i~k — Tth < 0, crosses), or whether the surface energy terms are 
omitted (triangles). 

feet of other forces i.e., \2(E t h + — Tth ~ r fc) + E m + r m )0 it is necessary 
to evaluate each term in addition to the sign of the surface terms. Since the 
evaluation of one type of energy is rather independent, it is hardly possible to 
obtain from a single observational technique all the relevant quantities necessary 
to perform this comparison. Thus, in a first approximation, one might attempt 
to measure the individual energy ratios (thermal/gravity, magnetic/gravity, and 
kinetic/gravity) separately. This information is available to us in the simulations 
and can be used to evaluate the corrections that should be applied to the ob- 
served classical indicators in order to take into account the effect of the surface 
energy terms (which represent the effect of the environment). 

2. Simulations 

In this study, we have analyzed a sample of 3D, isothermal, magnetized, driven, 
and self-gravitating MC simulations. The simulations vary by the strength of 
the initial magnetic field in the parent cloud. The MHD equations are solved 



1 E t h = | J v PdV is the volume thermal energy, where P is the thermal pressure, Ek = 
| J rriiV^dV the volume kinetic energy, E m = ^ J B 2 dV the volume magnetic energy, 
Tth = 5 f s TiPhidS the surface thermal energy, T k = ^ § s ripViVjhjdS the surface kinetic 
energy, r m = <j> s riTijhjdS the surface magnetic energy, where Ttj is the Maxwell stress tensor 

given by T = (1/4 tv) (B B — 0.5 B 2 T) and / the identity tensor, r is the position of each 
particle/cell with respect to the center of mass (CM), and v, the relative velocity of parti- 
cles/cells to the CM. The gravitational term is defined as W = — J pri{d<f)/dri)dV , where cj> 
is the gravitational potential. 



Figure 2. A schematic view on how to measure the kinetic surface energy 
term for a dense molecular core using a non-depleted high density tracer such 
as N 2 H + . A second lower density tracer (e.g., a CO line) is needed to obtain 
information about the gas dynamics outside the N 2 H + core boundary. The 
surface kinetic energy term is obtained by integrating the velocity gradients 
on the boundary. 

using the TVD code described in Kim et al. (1999) at a numerical resolution of 
256 3 . Periodic boundary conditions are used in the three directions. In physical 
units, the simulation box represents a scale of 4 pc, an average density of n = 500 
cm -3 , a thermal sound speed of 0.2 km s _1 , and turbulence is driven such as to 
maintain a sonic Mach number of 10 (i.e., the rms velocity is 2 km s" 1 ). The box 
has an initial Jeans number of 4. Results are presented at epochs that precede 
the appearance of runaway gravitational collapse in any core in order to avoid 
violating the Truelove criterion (Truelove et al. 1997). CCs are identified in the 
simulations by using a clump-finding algorithm based on a density threshold and 
a friend-of-friend criteria. The algorithm is described in more detail in Dib et 
al. (2006a). The sample of objects selected in each simulation is a mix of CCs 
identified at 5 density thresholds (7.5, 15, 30, 60, and 100 n). 

3. Results and Conclusions 

In the following, we show, as an example, the results for a model where the 
parent cloud in strongly supercritical, i.e., with an initial weak magnetic field of 
4.6 jiG and a mass-to magnetic flux ratio, normalized to the critical value for 
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gravitational collapse, of ^b ox = 8.8). In this model, we observe the existence 
of both collapsed, and other, non-collapsing, long-lived condensations (see Dib 
et al. 2006a for details). Fig. Q] (left) displays the relationship between the CCs 
ratios of thermal to gravitational energies (E t h — Tth)/\W\ and their respective 
(1/J^) estimates at the selected timestep. The two quantities are well correlated 
(stars in Fig. [[]). A linear fit is over-plotted to the data and is seen to follow: 



This relation provides the correction that ought to be applied to estimates 
of J c in CCs in order to get a more accurate estimate of the true thermal to 
gravitational energies ratio which takes into account the effect of the surface 
thermal energy term. If the surface terms are omitted (rth = 0, triangles in 
Figs. [1]), the linear fit shows, as expected, a slope close to 1 (0.99 ± 0.05). 
The comparison of the kinetic+thermal to gravitational energies ratios to a v i r is 
complicated by the fact that the quantity {E^+E^—Tk—r^ can be of either sign. 
However, for both positive and negative families, Fig. Q] (right) do not show any 
correlation between a v i r and the (kinetic+thermal) / gravity energy ratios (If the 
surface terms are omitted the slope is close to 1, i.e., 0.90 ±0.09, triangles). This 
result is essentially due to the complexity of the surface kinetic energy term, r& , 
which varies substantially from one object to another even when the CCs have 
a similar volume kinetic energy. Therefore, it seems likely that the evaluation 
of Tfc should be made directly from the observations. A simple method might be 
to integrate the kinetic energy for pixels which lie on or close to the boundary 
as defined arbitrarily by the molecular tracer and compare this quantity to the 
volume kinetic energy. A second method (Fig. [2D might rely on measuring the 
velocity gradients across the boundary of the CC. Such method requires high 
spatial resolution in order to sample the boundary and the use of two distinct 
molecular tracers such that the line of sight velocity measurements at the inner 
side of the boundary would not be affected by material which lies outside the 
boundary. Making some additional assumptions such as the isotropy of turbulent 
motions in the CC, it might be possible to use the velocity gradients in order 
to calculate Tt (the density at the boundary is equal to the critical density of 
the high density tracer). A detailed study, at higher numerical resolution, and 
including the comparison of the magnetic terms will be presented in a future 
work (Dib et al. 2006b in preparation) 
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